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Abstract 

Let d = (di, d 2 , ■ ■ ■ , d n ) be a vector of non-negative integers with even sum. 
We prove some basic facts about the structure of a random graph with degree 
sequence d, including the probability of a given subgraph or induced subgraph. 

Although there are many results of this kind, they are restricted to the sparse 
case with only a few exceptions. Our focus is instead on the case where the average 
degree is approximately a constant fraction of n. 

Our approach is the multidimensional saddle-point method. This extends the 
enumerative work of McKay and Wormald (1990) and is analogous to the theory 
developed for bipartite graphs by Greenhill and McKay (2009). 

1 Introduction 

Let d = (d 1 , d 2 , ■ ■ ■ , d n ) be a vector of non-negative integers with even sum. Let X = (xj k ) 
be a symmetric n x n matrix over {0, 1} with zero diagonal. Define G(d,X) to be the 
number ofnxn symmetric matrices A = {aj k ) over {0, 1} with zero diagonal, such that 

(i) row j sums to dj, for 1 < j < n; 

(ii) dj k = whenever Xj k — 1, for 1 < j, k < n. 

Equivalently, G(d, X) is the number of labelled simple graphs with n vertices of degree 
di, d 2 , ■ ■ ■ , d n , having no edges in common with the simple graph X. The special case 



'Research supported by the Australian Research Council. 



1 



where X is the zero matrix will also be denoted G(d). Define 
where a;,- is the sum of the jth row of X. 

One motive for interest in G(d, X) is that the ratio G(d, X)/G(d) is the probability 
that a random simple graph with degree sequence d has no edge in common with X. 
Similarly, G(d — x, X)/G(d) is the probability that X appears as a subgraph. In these 
cases, and throughout the paper, probability spaces have the uniform distribution. 

Define the matrix X = (xj k ) over {0, 1} with Xj k = 1 iff j ^ k and Xj k = 0. For 
convenience we will adopt the convention that Y^ijkex means the sum over all {j, k} such 
that Xj k = 1, and similarly Yljkex means the sum over all {j, k} such that Xj k = 1. Note 
that the equal sets {j, k} and {k,j} do not appear as separate terms in these sums. 

Define the following key parameters. 

(the number of edges) 

(the average degree) 
(the density ignoring the diagonal) 

(the number of edges of X) 



Direct asymptotic estimation of G(d, X) for nonzero X has been previously restricted 
to the sparse range. For representative results with bounded or very slowly growing 
degrees, see Bollobas and McKay [3] and Wormald [21]. For somewhat higher degrees we 
have the following. Let <i max = max^ dj, x max = max^ Xj and A = d max (d mAX + x max ). 

Theorem 1 ([15]). Suppose <i max > 1 and A = o(E). Then, as n — >■ oo, 

( ' } = exp v m — 

- E ^ xdA +0(A 2 /E) 
2E 

The error term in Theorem[T]is o(l) only under the stronger condition that A 2 = o(E), 
which implies that the graphs are quite sparse. The special case G(d) was determined by 
McKay and Wormald [20J under the weaker condition <i^ ax = o(E). 
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The probability G(d, X)/G(d) of being edge-disjoint from X, and the probability 
G(d — x, X)/G(d) of containing X as a subgraph are easily deduced from Theorem [TJ 
They can also be found directly over a sometimes wider range of d values. Let (a) b 
denote the falling factorial. The following is a consequence of Theorems 2.9 and 2.10 of 
McKay [1] 



Theorem 2 ([H]). If A + X = o(E) then, as n oo ; 

G{d- X ,x) _ n;=i(^- 



G(d) 2 X (E) X 



exp(0(AX/E)). 



In the case of dense matrices, the first asymptotically precise enumeration result was 
that of McKay and Wormald [19] , who proved Theorem [3] (below) in the case of X = 
with a slightly weaker error term. This has been extended by Barvinok and Hartigan pQ . 
They identify a symmetric matrix (A jfc ) over {0, 1} introduced in [19] (and used in greater 
generality in Section [3]) as the matrix that maximises a certain entropy function. They 
then express the asymptotic value of G(d) as an effectively computable function of (A jfc ) 
provided the values {Xjk}j^k are uniformly bounded away from and 1. This forces the 
average degree to be Q(n) but allows a much greater variation of degrees than we allow. 
They also show that (A,-^) matches a typical graph with degree sequence d in a sense that 
we will describe in Section [2J This theme is explored in a somewhat different way by 
Chatterjee, Diaconis and Sly [5]. 

Despite the absence of precise enumerative results for densities between o(n -1 ^ 2 ) and 
c/ logn, Krivelevich, Sudakov, Vu and Wormald [12] determined several almost-sure prop- 
erties of random regular graphs over various ranges of density using a combination of 
switchings and analysis. Other such properties were determined by Boldi and Vigna [2J, 
and Cooper, Frieze, Reed and Riordan [SIE]- 

More recently, Krivelevich, Sudakov and Wormald [UJ determined the probability 
of small induced subgraphs in random regular graphs of degree (n — l)/2 under some 
conditions on the order and degree sequence of the subgraph. 

The corresponding problems for bipartite graphs and digraphs were studied by Green- 
hill and McKay [8]; see that paper for a bibliography. The proof method in [8] is quite 
similar to that here. 

We will also have need for the following additional parameters, for 1 < j < n and 
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m>l. 



k\jkex 3=1 



R e. = <*j, -X* - X] ' 

j=i i=i 

jkex jkex 

n 

jk€X j=l 

K = ^(dj-dXd.-d). 

jkGX 

To calibrate and motivate our main enumeration result, we first develop a naive esti- 
mate of G(d,X) by extending an idea introduced in [19J. Generate a random graph by 
independently creating an edge jk with probability A for each jk G X. Each graph with 

E ( 71 1 X E 

E edges (none in common with X) appears with probability A (1 — A)^ 2 ' _ . Moreover, 
the event Ej that vertex j has degree dj has probability ( n ~^~ XjS j\ dj (1 — X) n ~ 1 ~ x i~ d 3 for 
each j. If we (incorrectly) assume that the events E\, . . . , E n are independent, we obtain 
a guess for G(d, X) as follows: 

G(d,X) = (1 - A)-*(A*(1 - A)^)(») { d ~ Xj )- 

In [19] it was proved that 



G(d) = V2G{d,Q) exp - 



R 2 



\QA 2 n 



+ o(i; 



under certain conditions on d. Our main result extends this to nonzero X. 

Theorem 3. Let a,b > be constants such that a + b < |. Then there is a constant 
e = e(a, b) > such that the following holds. Suppose that dj—d, Xj = 0{n l ^ 2+£ ) uniformly 
for 1 < j < n, that X = 0(n 1+2e ), and that 

n 

mm{d, n — d — 1} > 



3a log n 

for sufficiently large n. Then, as n — )■ oo ; 



G(d, X) = V2 G(d, X) exp ^ ^ 



0(n ) 
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Proof. The proof of this theorem is the main task of the paper. Here we will summarize 
the main phases and draw their conclusions together. The basic idea is to identify G(d, X) 
as a coefficient in a multivariable generating function and to extract that coefficient using 
the saddle-point method. In Section [3j we write G(d,X) = P(d, X)I(d, X), where 
P(d, X) is a rational expression and I(d, X) is an integral in n complex dimensions. 
Both depend on the location of the saddle point, which is the solution of some nonlinear 
equations. Those equations are solved in Section and this leads to the value of 
P(d,X) in (126]) . In Sections 13. 3113. 5] the integral I(d,X) is estimated in a superset of 
a small region 1Z enclosing the origin (defined in ( 128]) ) and an equivalent small region 1Z' 
enclosing (jr, . . . , tt). The result is given by Lemma [12] Finally, in Section 13.61 we note 
that the integral restricted to the exterior of 71 U 7Z' is negligible. The present theorem 
thus follows from (jl]) , (126]) and Lemmas [12] and HH 

The proof in various places requires e to be sufficiently small, but there are only a finite 
number of such places so we can choose e(a, b) to satisfy all of them at once. Note that 
the theorem remains true if e(a, b) is decreased, since the conditions become stronger. □ 

Throughout the paper, the asymptotic notation 0(f(n)) refers to the passage of n 
to oo. We also use a modified notation 0(f(n)), which is to be taken as a shorthand for any 
expression of the form 0(f(n)n C£ ) with c a numerical constant (perhaps a different con- 
stant at each occurrence). Under the assumptions of Theorem[3] we have A -1 , (1— A) -1 = 
O(logn). This implies, if ^,02,03,04 are constants, that A Cl (l — \) C2 n c ' i+C4£ = 0(n C3 ). 

2 Subgraph probabilities 

Define functions miss(d, X) and hit(<Z, X) as follows. The probability that a random 
simple graph with degrees d has no edges in common with X is 

(1 - A) x miss(d,X), 

and the probability that it includes X as a subgraph is 

X x hit(d,X). 

In this section, we apply Theorem |3] to estimate these probabilities. To avoid unnecessary 
messiness regarding the value of s, we will suppose e(a, b) in Theorem |3] is chosen to be 
small enough to satisfy the finite number of places in this section where a statement is 
only true if e is small enough. 
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Theorem 4. Under the conditions of Theorem^, we have 



XX XX 2 X(1-2X)X 3 XX 2 D 



miss(d, X) = exp — — — — h — — — — — I 2 2 + 



(l-A)n 2(1-X)n 6(1-A)V (l-A)n 2 A(l-A)n 2 
Ci,i (1— 2A)C 12 C 2 i _ 6 , 



(1 - A)n 2(1-A)V 2(1-A)V 



and 



Mttd.X) = exp| S=^L - (1±^2 _ (Wffg. + (1^ 



An 2An 6An An 

+ n< 2 2 TTnm + 0(n ) . 



A(l-A)n 2 An 2A 2 n 2 2A 2 n 2 

Proof. Since (1 — X) x miss(d, X) = G(d, X)/G(d), the first part can be obtained from 
Theorem [31 The second part can be found in similar fashion, or by noting that the 
probability of a random graph avoiding X is the probability of the complement of the 
graph having X as a subgraph. □ 

If X is not too dense, the probabilities in Theorem H] asymptotically match those 
for an ordinary random graph with edge probability A. Sufficient conditions are that 
miss(d, X) = 1 + o(l) if 

AX 2 + Xmax, |<L — d\ = o((l — A)n), 

and hit(d,X) = 1 + o(l) if 

(1 — A)X 2 + X maXj \dj — d\ = o(An) . 

Both these sufficient conditions hold, for example, if X = O^n 1 ^ 2 ^ 26 ), or if dj—d and Xa 
are uniformly 0{n e ) for 1 < j < n and X = 0(n 1_2e ). 

Since Theorem 0] is rather complex, we give some special cases to facilitate its applica- 
tion. We also give the value of num(d, X), which is the exponential factor in Theorem El 

Corollary 5. Suppose the conditions of Theorem [3] hold, and in addition assume that 
d 1 = ■ ■ ■ = d n = d. Then 

1 X(X 2 -H) 
4 + (l-A)n 2 
AX AX 2 A(2-A)X 3 



num(d, X) = exp ( - H — ^ — - - 2 ; + 0{n 1 
miss(d, X) = exp 



x 2 2 
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hit(d, X) = exp 



;i_ A )x (1-X)X 2 (1-X 2 )X 3 



Xn 2An 6A 2 n 2 

+ (!z^_(i = J0£ + o( „- u 



An An 

Corollary 6. Suppose the conditions of Theorem [3] hold, and in addition assume that 
%i = ■ ■ ■ = x n = x (which implies that x = 0{n 2e ) ). Then 



i -\r i ' K R „- _/, 

^ {d ,X) = e^[- + w -^-— 2 -^-, + 0(n 



/ Xx(x-2) xR K _ b , 

m 1S8 (d, X) = exp( - - — 2 + 0(n 

, / (1— A)x(x— 2) xR K ^. _ b . 

h,t( d , X) = exp( ' 4A ' - ^ - ^ + 0(» 

The two parts of Theorem H] have a common generalization. Let Y be a supergraph 
of X. Then the probability that a random graph with degrees d has intersection with Y 
equal to X is 

G(d-x,Y) 



G(d) 

If the degrees of Y are y x , . . . , y n , with = 0(n 1 ^ 2+e ) uniformly over j, and 527=1 Vj 



0(n 1+2e ), then this probability can be computed using two applications of Theorem [3j 
The resulting general formula is rather complex, so we will be content with presenting 
the special case where Y consists of a single clique and otherwise isolated vertices. This 
is the important case of an induced subgraph. 

Suppose that for some m, we have x m+ i — ■ • • — x n — 0. Let X™ be the subgraph of 
X induced by vertices 1, . . . , m (so x' m ' has the same edges as X). For k, £ > 0, define 
the quantity 

m 

u k ,e = ^2( d j ~ d) k (xj - X(m - 
i=i 

Theorem 7. Assume the conditions of Theorem^ and in addition that m = 0{n l ^ 2+£ ) 
and x m+ i = ■ ■ ■ = x n = 0. Then the probability that a random graph with degree sequence d 
has x' m ' as an induced subgraph is 



A (1 — A)^ : 

(2u lx -u Q2 m (l-2A)a; 01 4^ a; 0il - a; 2 ^ - 2w 2 

x exp — H 1 H ■ ■ 7j 

\ AAn 2n AAn 8An 

(2^! x - u 2 o - u )2 )m (l-2A)(w 3 + 3w 21 -3w 12 ) , „, _ 6 s 

H ■ 7-5 ■ ; 5-V ; VOin ) 

AAn 2AA 2 n 
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Note that, within the stated error term, the probability is independent of d m+1 , ...,d 
except inasmuch as they contribute to d and A. 

The factor A (1 — A)^ 2 ' in Theorem [7J is the probability for an ordinary Erdos- 
Renyi random graph with edge probability A. A sufficient condition for the argument of 
the exponential to be o(l) is m 2 (m + max™^^ — d\) = o(An). Relaxing this condition 
by a factor of m allows us to see the leading terms of the deviation of behaviour from an 
ordinary random graph. 

Corollary 8. Assume the conditions of Theorem^ and also thatm(m + max.™ =1 \dj — d\) = 
o(An) and x m+1 = ■ ■ ■ = x n = 0. Then the probability that a random graph with degree 
sequence d has x' m ' as an induced subgraph is 

\Xf-x wM-X I W 0,2 . 

(1 " A)UJ eXP ^ _ ^ + 0(1) j- 

The term a; 0)2 / (4An) was obtained in [11] in the case that d = ((n— 1)/2, . . . , (n— 1)/2), 
m = o(y/n), provided x 1} . . . ,x m don't differ too much from Am. (In the regular case 
u 1A = 0.) 

It has been shown by Barvinok and Hartigan [1] , see also [5] , that an independent-edge 
model more accurately matching the d model has each edge jk chosen with probability Xj k , 
where these constants were introduced in [19] and will appear generalised in the following 
section. We will call this the {Xj k } -model. Under our strict constraints on d, we have 

a^-d 4-d (l-2X)( dj -d)(d k -d) ~ / 
1 n n 2An 

We can restate Theorem [TJ with that model in mind. 

Corollary 9. Assume the conditions of Theorem^ and in addition that m = 0{n l ^ 2+e ) 
andx m+1 — ■ ■ ■ — x n — 0. Then the probability that a random graph with degree sequence d 
has X™ as an induced subgraph is 

n ^ n c 1 - ^ 

jk&X jk<£X 

^0,2 , m (l-2A)w i , 4Ux nW ,1 - ^o,i 



x exp ( + — H — - + , 2 

, — , -—1,2) — i)\ 

+ nA ,2 2 + ) 



4An 2n 4An 8An 2 

(2a; M - u 0)2 )m (1-2A) (w„ i3 - 3a; 1)2 ) 



4An 2 24AV 



where the two products are restricted tol<j<k<m. 
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Note that the uj x ^ term of Corollary [8] has disappeared, but the w 2 term remains. 
The exponential factor quantifies how much the {A jfc }-model is in error. However, when 
X is generated according to the {A jfc }-model, the expectation of the argument of the 
exponential is 0(n~ b ) and the variance is Ofa^ 1 ^), so there is some sense in which we 
can say that the {A jfc }-model gives a very accurate estimate for typical subgraphs. The 
details of this remain to be worked out. 

Let Y be any graph on n vertices and let Y be its number of edges. Define the random 
variable E dY to be the number of edges that a random graph with degree sequence d 
has in common with Y . Barvinok and Hartigan [lj proved that E dY is concentrated 
close to ^2j keY Ajfc when Y = £l(n 2 ). In fact their estimate is explicit enough to infer 
this concentration with weaker bounds for Y = u{rC'^ 2 logn). Under our strict conditions 
on d, we can obtain such a result for all Y. We have not determined the best result 
that follows from Theorem [3j but will for this paper be content with the following weak 
corollary of Theorem [7J 

Corollary 10. Assume the conditions of Theorem^ If Y = 0(n 1//5 ) ; then 

Prob(E dX ) = k) = ^A fc (l-A)(")- fe (l + 0(n-^)) ( 3 ) 

uniformly over k — 0, . . . , Y. If Y = fi(n 1//5 ) ; then 

E dX = \Y{l + O{n- 1 ' 10+s )) 
with probability 1 — 0(exp(n~ 5 )) for any 5 £ (0, 4). 

Proof. If Y = 0{n 1/5 ), there are 0(n 1/ ' 5 ) vertices incident with edges of Y, so Theorem [7J 
shows that Z has a binomial distribution with the precision given by (J3j). If Y = ^(n 1 ^ 5 ), 
divide the edges of Y into subsets of size 0(n 1//5 ) and apply (j3J) and Chernoff's Inequality 
to bound the number of edges in each subset that are in common with the random graph. 
(Clearly the constants in this corollary can be tuned in various ways.) □ 

The theorems above should be enough to allow transfer of quite a lot of the theory 
of ordinary random graphs to dense random graphs with given degrees. However, our 
purpose in this paper is to develop the tools rather than to explore the applications in 
detail. We will be content with some simple illustrations. 

Theorem 11. Let d = (d, d, . . . , d) satisfy the conditions of Theorem^ Then for a ran- 
dom d-regular graph, we have the following. (Note that in each case the quantity in front 
of the exponential is the expectation for ordinary random graphs with edge probability X.) 
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(a) If n is even, the expected number of perfect matchings is 

A" /2 n! /1-A _ 6> 
exp I — — — h (J[n , 



2 n/2 (n/2)\ \4A 

(b) If q = q(n) is a integer function such that 3 < q < n, then the expected number of 
q-cycles is 

avj ffip f_(hiM +0( , r , : 



2q(n—q)\ \ Xn 
(c) The expected number of spanning trees is 

n-2-vn-l ( — A) _ 6 , 

n A exp I — — V U[n 

\ 2 A 

Proof. Parts (a) and (b) follow immediately from Corollaries |6] and [5j respectively. 

Part (c) is not so simple since trees have various degree sequences and those with 
maximum degree greater than n 1//2+e do not satisfy the requirements of Theorem [31 Let 
T be the set of all labelled trees with n vertices. If x denotes the degree sequence of a 
member of T, let 7[ be the subset of T with x max < n £ and let T 2 —T\7t- 

For sufficiently small e > 0, the following are true. 

(i) The probability in T that the maximum degree exceeds k is at most 2n/k\ for any 
integer k > 0. 

(ii) X = n-1. 

(iii) In T\ we have that X 2 = 5n + 0(n 1//2+3e ) with probability 1 — 0(e n ). 

(iv) In 7i we have X 3 = 0(n) and if = 0(n). 

(v) The sum of A -E" =1 max{o^-n E } Qver ^ ig ^ n »-i/( n£ )i _ 

Facts (i) and (iv) follow from the well-known generating function for labelled trees by 
degree sequence, which is 

z x z 2 ■ ■ ■ z n {z x + z 2 + • • • + z n ) n ~ 2 . 

To obtain (iii), note that the same probabilities occur if we take Xj = 1 + Yj for 
1 < j < n, where Y 1} . . . , Y n are independent Poisson variates with mean 1 truncated at 
n , subject to having sum n — 2. (This is a standard property of multinomial distributions; 
any mean will do.) Now we can write 



Probf \YXf -2n\> n 1/2+3s YXj =n-2)< 



Probfl^r/ -2n\ > 



l/2+3e 



Piob(j2Y j = n-2\ 
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The expectation of Y 2 is 2 + 0(rT l ). Now bound the numerator by applying a concen- 
tration inequality like Hoefffing's [10] and the denominator by noting that Yj would be 
a Poisson distribution with mean n except for the truncation. Item (iii) follows. 

To obtain (v), note that YTj=i m ax{0, Xj — n £ } > A — n e for trees with maximum 
degree A, and the number of such trees is bounded by (i) with k = A — 1. Summing over 
A > n £ gives the desired bound. 

Now we bound the expected number of trees in T 2 that appear in a random <i-regular 
graph. For a tree T e %, let F(T) be any forest obtained by deleting all but \n e ] edges 
from each vertex that has degree greater than n £ . Then the probability that T appears is 
bounded by the probability that F(T) appears. Moreover, F(T) satisfies the requirements 
of Corollary [5] and has at least n — 1 — YTj=i max{0, Xj — n e } edges. Corollary [5] gives 
hit(d,F(T)) = e° (nE) . Applying fact (v), we find that the expected number of these trees 
easily falls within the error term of part (c) of the theorem. 

Finally, the trees in 7i all satisfy the conditions of Corollary |5] and have hit(e£, X) = 
0(n 3 ). Those with X 2 = 5n + 0(n 1/2+3e ) have hit(d,X) = exp(7(l-A)/(2A) + 0(n" b )). 
Part (c) of the theorem now follows from items (i) and (iii). □ 

The average number of spanning trees in random regular graphs of bounded degree 
was studied in [13]. 



3 Proof of Theorem 3 



In this section we express G(d, X) as a contour integral in n- dimensional complex space, 
then estimate its value using the saddle-point method. 

We will use a shorthand notation for summation over doubly subscripted variables. 
From the matrix X = (xj k ), define sets 

X{j) = { k : 1 < k < n, x jk = 1 }, X(j) = { k : 1 < A; < n, x jk = 0, k ± j } 

for 1 < j < n. Note that j X(j), X(j), also that |X(j)| = Xj and |X(j)| = n — 1 — Xj. 
If Zjk is a symmetric variable for 1 < j, k < n, we define 



3* 



/] z jki z ** z jki 

k=l j=l k=l 
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kex(j) jkex 

z j*\x = z i k) z **\ x = z i k ' 

kex(j) jkex 

There is some slight lack of symmetry in the definitions. To clarify, we note that 
^2 z j*\x = 2z**\x and ^2 z j*\X = 2z **\Xi but 



3 



Firstly, notice that G(d, X) is the coefficient of z 1 1 z 2 2 ■ ■ ■ z n n in the function 



Y[ (1 + ZjZ k ). 

jkex 



By Cauchy's theorem this equals 



1 / / Ujkex (! + z j z k) 



where each integral is along a simple closed contour enclosing the origin anticlockwise. 
It will suffice to take each contour to be a circle; specifically, we will write 

Zj = Tj e > 

for 1 < j < n. Also define 

1 + rf k 

for 1 < j, k < n. Then 

G(d,X) = n ^ x _ \ k> ... n ^ x{ ^ ^dO, (4) 



where = (9 1 ,...,d n ). Write G(d, X) = P(d,X)I(d, X) where P(d,X) denotes the 
factor in front of the integral in (jlj) and I(d, X) denotes the integral. We will choose the 
radii Tj so that there is no linear term in the logarithm of the integrand of I(d, X) when 
expanded for small 0. The linear term is 

n 

Yl x jk(°j + e k) -z^dj^j. 

jkeX 3=1 

For this to vanish for all 6, we require 

\ j .\X = d j (l<j<n). (5) 
12 



Although it is not hard to show that ((SJ) has an exact solution, we can get by with a near- 
solution since (]1J) is valid for all positive radii. In Section 1331 we find such a near-solution 
and determine to sufficient accuracy the various functions of the radii, such as P(d,X), 
that we require. In Section [3T31 we evaluate the integral I(d, X) within a certain region 1Z 
defined in (128]) . Section [3T61 notes that the contribution to the integral from the region 
outside of TZ and its translate 1Z + (tt, . . . , it) is minor in comparison. 



3.1 Locating the saddle-point 

In this section we derive a near-solution of (jSD and record some of the consequences. As 
with the whole paper, we work under the assumptions of Theorem [3j 
Change variables to {aj}j =1 as follows: 

r j = r / + 2 J » ( 6 ) 
1 — r a,j 

where 



A 



From ([6]) we find that 



where 



1-A 

X jk / X = l + a j + a k + Z jk , (7) 



ctja k (l — r — r dj — r a k ) 

Z jk = " ' 2 ' 

1 + r ajd k 

and that equation (jSJ) can be rewritten as 



-j- = (n- l)aj - a s Xj + a k + Z j*\x- ( 9 ) 



A 

kex(j) 

Summing ([9]) over all j, we find that 

n 

X = ~ a i x i) +Z **\x- ( 10 ) 

i=i 

Replace the term ^2 k& x(j) a k m & by X/fc=i a k~ J2kex(j) a k~~ a ji an d substitute the value 

n 1 n XI 

/ , a fc = — / , ( a fc + a k x k ) H 

fe=i fe=i 
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implied by (fit)]) . After some rearrangement, we find that aj = Aj(ai, . . . , a n ) for each j, 
where 

A,(a x , . . . , a n ) = 1 2 2 / ^l a fc + ^k-^k) 

An n n n ,^ 

1 1 1 ( U ) 

+ ^ E ^-- Z i*|x + ^HX- 

In the vicinity of a = (0, 0, . . . , 0), the iteration a := (Ai(a), . . . , A n (a)) is a contraction 
mapping that converges to a solution of as can be proved using the method demon- 
strated in jl]. However, as noted above, we do not need to solve ([5]) exactly but will 
work with an approximate solution. Hopefully without confusing the reader, from now 
on we will use a to denote the result of four iterations starting at a = (0, 0, . . . , 0). We 
will also write Zj k and Xj k to mean the values implied by ([7]) and OH]) for our chosen a. 
Applying ( ITTj) four times, we find 

' = ^ + r^-^ + r^ + '" +<5( ' r5/2) - (12) 

au \n n An 

where the ellipsis conceals about 60 terms of order 0(n~ 3//2 ). Most of the terms involve 
counts of subgraphs of X up to order 5, with the vertices weighted by powers of the 
numbers {Sj}. This implies an expansion 

3 2 A An V ; 



The value of Aj k is given by substituting estimate ffT21) into ([7]). In particular, uniformly 
over j, 

VP 

Define a jk , jk , >y jk by a jk = jk = >y jk = if j = k and 



\ H x=d J + 0(n- 3 / 2 ). (13) 



~~ ^jk) — A + a jk , 
±\ jk (l - X jk )(l - 2X jk ) = A 3 + p jk , (14) 
— X jk )(l — 6\ jk + 6A 2 fc ) = A A + 7 jfc , 



if j ^ k, where 



A = ±A(1 - A), A 3 = lA(l - A)(l - 2A), and A 4 = ^A(l - A)(l - 6A + 6A 2 ). 
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In evaluating the integral I(d,X), the following approximations of aj k , f3j k , and jj k 
for j ^ k will be required: 



(l-2A)((5 J +4) 5 2 + 5 2 k (l-12A)fy^ (l-2A)(^ + ^ fc ) 
2n 2n 2 4An 2 2n 2 

_ A(1-2A)X (l-2A)(^ + ^x fc ) 5(ri -3/ 2) 
n 2n 



3.2 Estimating the factor P(d,X) 

Let 

A = n ^(i-a,,) 1 - 

jfcex 



Then 



jkex 

n 



n ( i+r i r fc) n r / 3 *' 

jk&X i= l 



n c 1 + r ^) n 







i+0(n" 3/2 ) 



(15) 



(l-lM)(tf J + fr) + ~ (16) 



Pjfe ,, 
on 

Tife = 0(n" 1/2 ). (17) 
We will also need the following summations. 

a, = 1(1 - 2A),, - | - * + (1 - 2A '^f + B ' ) + 5(„-V») (18) 

a,. = -— + A(l - 2A)X + 5(n 1/2 ) (19) 
n 

^ = ±(1-12^ + 6(1) (20) 

= O(n) (21) 



using (|T3l) . Therefore the factor P(d, X) in front of the integral in (Tj0) is given by 

P(d,X) = (2ny n A- 1 exp(d(n^ 1/2 )). (22) 
15 



We proceed to estimate A. Writing Xj k = A(l + Zj k ), we have 



A A (1-A)^ A ) ~^^\\-X 



(23) 

+ 7m T\ z jk ~~ 777; WI Z % "I V , TT3 — ~ z jk + 0(n 5//2 ). 



A 2 A(l -2A) 3 A(1-3A + 3A 2 ) 
2(1- A) S ' fc ~ 6(1 - A) 2 Zjk + 12(1 - A) 3 * ik 



We know from (TT3~j) that A^jf = £' + 0(ra~ 1 ^ 2 ), which implies that = X + 0{n^ 1 ^ 2 ), 

hence the first term on the right side of (T231 contributes A AX (1 — A)~ AX exp(0(n -1 ^ 2 )) 
to A. Now using ([7]), and recalling that |X| = (™) — X, we can write Zj k = ctj + a k + Zj k 
and apply the estimates in the previous subsection to obtain 



A 



(A A (1-A)^ A )^(1-A)" 



x 



(n + 2)R 2 (1-2X)R 3 (1-6A)R A 
X6XP( AAn 2 24AV + 96A 8 n 8 ( 24 ) 



C 21 + 2D A 2 X 2 i? 2 ~ _ 1/2 

9 9 H O-^+^l^ ) 

4An 2 2An WA 2 n V 



As in [19], our answer will be simpler when written in terms of binomial coefficients. 
Using Stirling's formula or otherwise we find that 



-j-j- fn Xj 1\ = ^ n yn/2^-n/2-Xn(n-X)Q _ ^-n/2-(l-A)n(n-l)H 



-2X 



I 1-144 R 2 (1-2X) 2 R 2 (1-2X)R 3 , nK , 

x exp — — - + ^ + v ^A^- 25 

• { 2AA AAn 16A 2 n 2 2AA 2 n 2 

_ (1-6A)R A C 2tl XX ~ 1/2 , 

96AV AAn 2 (l-A)n 1 

Combining (|22]h (plj) and (j25j) ; we find that 

-n/2 



l-A 



i=i 

x exp 



1 - UA R 2 R\ 



2AA 16An 16AV 

AX XX 2 D _ 6 , 

(l-A)n + (l-A)n 2 2An 2 + ^ ' 



(26) 
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3.3 Estimating the main part of the integral 



Our next task is to evaluate the main part of the integral I(d, X) given by 

I(d,X)= ■■■ ll ^ x{ >ld6. (27) 

J-n J-* exp(z 2^ j=1 djVj) 

It will be established in this section and the next that the value of the integral is 
concentrated near the places where the integrand has the largest absolute value. This 
happens at the two points = (0, 0, . . . , 0) and 6 = [tt, tt, . . . , 7r). These two points are 
equivalent, since the integrand is unchanged under the mapping i->- + (tt, 7r, . . . , 7r). 
(This requires the fact that dj is even; otherwise the mapping changes the sign of the 
integrand and the integral is zero as it should be.) Consequently, in this section we will 
focus on a neighbourhood of (0, 0, . . . , 0), specifically the hypercube 71 defined by 

11= {0 : 1^1 <n~ 1/2+£ ,l <j <n}. (28) 

Let F(0) be the integrand of ( 12 7p . We are going to establish the following. 

Lemma 12. Under the conditions of Theorem^ there is a region S with 41Z 
such that 

n/2 



(29) 



Js v ; \AnJ F V 2AA (l-A)n lQA 2 n v 

In a region 0(1)7?., we can expand 

jkex jkex 

+ J2( A * + 7i*)(*i + + d(n- 1 ' 2 + \°j + h \ 

jkex jkex 

= exp(- ( A + a jk )(6 j + 9 k ) 2 -i ( A 3 + Pjk)(0j + e . 

l<J<fc<n l<j<fc<n 

+ Yl M 9 i + ^) 4 +Y A Pi + + 0(n- 1/2 )\ , 

l<j<k<n jkeX ' 

where A, A 3 , A 4 , ctj k , (3j k , and jj k were defined in (fHjl . Approximations for aj k , (3j k , jj k 
were given in (|T5|) -( 1T7|) . Note that otj k , (3j k , jj k = 0{n~ 1 ^ 2 ) uniformly over j, k. 

We will transform the integral to diagonalize the quadratic terms, proceeding in two 
steps. The first step will diagonalize the quadratic form Yli<j<k<n(®3 + ®k) 2 \ an d the 
second will complete the diagonalization. 



3 

k) 
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3.4 First change of variables 



We first adopt from [IH] a linear transformation that diagonalizes the quadratic form 
Ei<j<fc<n(0j + °kf- Define c and y = (y 1: y 2 , ...,y n )hy 



C =l-W^)=l-2- 1/2 + 0(.- 1 ) (30) 



<>, //, -T,H (!<i<")- (31) 



k=l 



The transformation = T^y) defined by f l3T|) has determinant 1 — c. Also 

(1 + c)K C Tf ^ C (1 - c)" 1 ^. (32) 
For £> 1, define //^ = E™ =1 2/j- We find the following translations. 

$^ = (1-0)01 

E (^ + ^) 2 = (^-2)/i 2 

l<j<fe<n 

E + ^) 3 = (n - 4)/i 3 + (3(1 - 2c) + 12c/n)/^ 2 
1 ^' < ^ n + ((-6c + 12c 2 - 4c 3 )/n - 4c 2 (3 - c)/n> 3 

E + ^) 4 = (n - 8)/i 4 + 3^ + (4(1 - 2c) + 32e/n)/^ 3 
1 - Kfc - n - (24c(l - c)/n + 48c 2 /n 2 )/i 2 /i 2 

+ (8c 2 (l - c)(3 - c)/n + 8c 3 (4 - c)/n 3 )/4 

E + e kf = E(( X ~ 4c / n ) a i* + 2c 2 a„/n 2 )?/ 2 

l<j'<fe<n j 

+ E ("j fc ~ 4c "i*/^ + %<?<x«*/n 2 )yjyk 

E + 9 hf = E( X " 6c / n + 12c V" 2 )/3i* - 4c 3 /3„/n 3 ) 2/ 3 



l<j<k<n 



+ E'((3 - Uc/n)P jk - 6c(l - 4c/n)/V™ 

j,k + I2c 2 (3 3 jn 2 - 12c 3 pjn 3 ) yj yt 

+ E (~ 6c /W™ + 12c 2 /3 j+ /n 2 - Ac z f3^/v?)y j y k y i 



3,k,i 



4c 4c 2 



E (fy + ^) 2 = E + ~ — ^ E + — • 
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In the above, and following, a summation is over 1,2, ... ,n for each index unless otherwise 
specified. Moreover, a prime on the summation symbol (as ^2 ) m eans that only terms 
where the summation indices have distinct values are included. For example, 

means 

j,k l<j'<n,l<A;<n 

Using the size of the hypercube 1Z together with the bounds (IT5|) - (|2"T|) . we find that 
whenever 6 G 0(1)7?., we have F(0) = G(y), where 



G(y) = - J2(( n ~ 2 ) A + a i* - Ax i)y 2 j 



2 

i 

2 



+ Axj k — 2Acxj/n — 2Acx k /n + AAc 2 X/n 2 )yjy k 
- i^(nA 3 + /3j*)yj 



- i J^(A 3 (-6c + 12c - 4c 3 ) jn - 6cf3 jk /n + 12cf3 j Jn) 

3,k,i 

+ nA 4 ^ Vj + 3A 4 Yl y2 j y l 

+ 4A 4 (1 - 2c) yj yl - 24A 4 c(l - c)/n y jVk y 2 

j,k j,k,£ 

+ 8c(l - c)(3 - c)A 4 /n ^ ViVkViVm + d{rT l/2 ). 



(33) 



j,k,£,m 



3.5 Completing the diagonalization 

We now make a second change of variables, y = T 2 {z), that diagonalizes the quadratic 
part of G(y), where z = (z l: . . . , z n ). We will use the method from [8J that is a slight 
extension of [TS] Lemma 3.2]. 

Lemma 13. Let U and Y be square matrices of the same order, such that U~ l exists 
and all the eigenvalues of U~ X Y are less than 1 in absolute value. Then 

(i + Yu^y 1/2 (u + y) (i + u^Y'y 1 ' 2 = u, 

where the fractional powers are defined by the binomial expansion. □ 
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If we also have that both U and Y are symmetric, then (J + YU 1 ) l ^ 2 is the 
transpose of (I + U~ 1 Y)~ 1 ^ 2 , as proved in ^j. 

Let V = (vj k ) be the symmetric matrix such that the quadratic terms of ( 1331) are 
—yVy T . We have for all j ^ k that 



Vjj = An + 0(n 1 ' 2 ), 
v jk = Ax jk + d(n~ 1/2 ). 

Apply Lemma [TBI with V = U + Y where U is the diagonal matrix with the same diagonal 
entries as V. The matrix U~ 1 Y has jfc-entry equal to n~ 1 Xj k + 0(n~ 3 ^ 2 ). Therefore, since 
the oo-norm (maximum row sum of absolute values) of U~ 1 Y is 0{n~ 1 ^ 2 ), the eigenvalues 
of U~ l Y are all d(n~ 1/2 ). 



Let T 2 be the transformation given by T 2 (y) = z, where z T = (J + U 1 Y) 



-1/2 T 

By jU Lemma 2], the Jacobian of T 2 is 1 + 0{rT 1 ^ 2 ). Expanding (I + U~ 1 Y)~ 1 ^ 2 we find 
that for y G 0(1)11, 

n 

Vj = z 3 + ^2{0(n~ 3/2 )z k + d(n~ x )xj k )z k , (34) 

k=l 

for each j, where the coefficients are uniform and independent of z. An expression of 
identical form writes z in terms of y. For y G 0(1)7?., we find that G(y) = H(z) where 

3 



)Z ;j 



H(z) = - ^2((n-2)A + a jt - Axj)z 2 - i^{nA z + 0(n 1/2 )), 
j j 

- i ^ (3A 3 (1 - 2c) + d{n~ xl2 ))zjz\ - i ^ 0[n~ x )zjZ k z l 
j,k j,k,£ 

+ nA 4 + 0(n -1 ))^ + 3A 4 ^'(l + 0(n' l ))z 2 z 2 k 

j j,k 



(35) 



+ 0(l)z jZ 3 k -J2' °( n 1 )zj*k4+ ^ 0(n 2 ) Zj z k z e z m + 0(n 1/2 ), 

j,k j,k,£ j,k,£,m 

with only the final expression of the form 0( ) being a function of z. 

Now define S = T 4 1 {T 2 1 (2TZ)). By (|32) and (I51jh K C S C 4K. Consequently the 
conditions for our approximations are satisfied and f[3"5"j) is valid for z G 2TZ. 

We can now apply Theorem [TH] (see Appendix) to estimate the integral of H(z) 
over 21Z. We list the coefficients required. 

A — A D jkt = 0(1) 

N = n Ej = A A + din' 1 ) 

Jj = F jk = 3A 4 + Oin- 1 ) 
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aj = 


(2A- 


- a J# + Axj)n 1 ^ 2 


Gjk 


= 0\ 


' -1/2 


Bj = 


-iA 3 


+ 0{n- 1 ' 2 ) 


Hjki 


= d[ 


' -1/2 


Cjk = 


-3(1 


- V2)iA 3 + 0{n- 1/2 ) 


I jklm 


= d( 





We can take A = | and have S(z) = 0{n l ^ 2 ). Applying Theorem [T5| we find that 

tt, n . f W 2 / 1 - 20A AX i? 2 ~ 1/2 . A \ . , 

#(*)<**=( — ) exp + + — JL- + 0( n 1/2 Z), 36 

where 

From the conditions of Theorem [3], we have that 0(n _1 ^ 2 )Z = O^ -6 ). Lemma [121 now 
follows on recalling that the Jacobian determinants of T l and T 2 are + 0(n -1 ^ 2 ) and 
1 + 0(n -1 ^ 2 ), respectively. 



3.6 Bounding the remainder of the integral 

In the previous section, we estimated the value of the integral I(d, X) restricted to a 
small region S D 1Z. As mentioned earlier, the integral over «S + (tt, . . . , 7r) is the same. 
It remains to bound the integral over the remaining parts of [—tt, vr] n . Define 1Z C = 
[—TT, IT (tt, . . . , 71")). By employing the same technique as in [8], but with the 

dissection of the region utilised in [19], we can establish the following. We will omit the 
proof since no new techniques are required, but we note for convenience that 

\ F (°)\ = II MA + a^Oj + Ok) 

where 

fjk(q, z) = v 71 - 4 ^! - cos z) < exp(-qz 2 + j^qz 4 ). 
Lemma 14. Under the conditions of Theorem^ 

[ \F(0)\dO = 0{n~ 1 ) [ F(0)dd. 
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Appendix: The value of an integral 



In this appendix we give the value of a certain multi-dimensional integral. A similar 
integral appeared in [19J and variations of it appeared in [9| [T6T4T8] . 

This appendix is notationally independent of the rest of the paper. Summations 
without explicit limits are over 1,2, ... ,N for each of the summation indices. A prime 
on the summation symbol (as Yl' ) indicates that only terms with distinct values of the 
summation indices are included. 

Theorem 15. Let e , e , e , e, e, A be constants such that < e < e" < e" , e > 0, s > 0, 

and < A < 1 . The following is true if e" and e are sufficiently small. 

^ / 

Let A = A(N) be a real-valued function such that A(N) = Q(N~ £ ). For 1 < j, k, £, m, 

let hj, Bj, Cj k , Dj ke , Ej, Fj k , Gj k , Hj ke , lj k i m , and Ja be complex-valued functions of N 

such that aj, Bj, . . . , J,- = 0(N £ ) uniformly over j, k, i, m. Suppose that 

f( Z ) = exp Lan e A + E fas + Nl/2 E ^ +N H + E' d dk *A 

3 3 3 3 j,k 

+ N' 1 Djkt z j z k z e + N E Ej z j + E Fjk z j z l + N 1/2 E Gjk z j z l 

3,k,£ j j,k j,k 

+ N Hjict ZjZ k Z£ + N ljkim z j z k z £ z m + $( z ) j 

j,k,£ j,k,£,m 

is integrable for z = (z±, z 2 , ■ ■ ■ , z N ) G U N and S(N) = max zeC / |5(;z)| = o(l), where 

U N = {z C R N : \zj\ < N~ 1/2+i for l<j< N}, 

where i = s(N) satisfies e < 2e < e" . Then, provided the OQ term in the following 
converges to zero, 

N/2 



p. / \JV/z 

/ f(z) dz =(-£- ) exp(<9 1 + 0{N- 1/2+s + (N- A + 6(N))Z)), 
Ju N \AN J 



where 



®i = — * — no °7 "i — ^ — "i — ^ — "i — ^ — o BAC jk 

2AN i/2A^ J AA N 3 16A 3 N^ 3 8A 3 N 2 ^ 3 3 

3 3 3 3,k 



+ 1QA 3 N 3 5, + ^ 2 N ? Ej + AA 2 N 2 ^ Fj 

3,k,t 3 3,k 

+ T77 E ^ + 7T^7 E + JITZZ E A'.fc^fc 



AiV^ J 4A Z N^ J J AA N 

3 3 3,k 
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Z = exp 




^ Im( ^' )2 + ^ Im( ^ Im( ^' fc) 

3 j,k 



-^^^ lm(B jk )lm(J k )J. 



Proof. The method of proof is the same as in [I], with extra terms added. To simplify the 
process, we did not explicitly compute the lower order terms which are presented as 2 
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